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3 ' Abstract 

. P . 

The operator associated to the angular part of the Dirac equation in the Kerr- 
Newman background metric is a block operator matrix with bounded diagonal and 
unbounded off-diagonal entries. The aim of this paper is to establish a variational 
principle for block operator matrices of this type and to derive thereof upper and 
<**j \ lower bounds for the angular operator mentioned above. In the last section, these 

Oh 1 analytic bounds are compared to numerical values from the literature. 

1 Introduction 

Variational principles are an important tool to give estimates for eigenvalues of a 
selfadjoint operator A since no knowledge about the corresponding eigenvectors is 
needed. The classical variational principle based on the Rayleigh functional ^p^r ; 
■ V> £ ^i-A) \ {0}, applies only to semibounded operators, see, e.g., [RS78]. 

For example, the eigenvalues of A below its essential spectrum are given by 

O ■ K = nun max . (1) 

, dim Li—ri 

o 

In this paper, however, we are interested in analytic bounds for the eigenvalues of the 
angular part of the Dirac operator in the Kerr-Newman metric which has been studied 
r> \ for instance in [BSW05] . It is neither bounded from below nor from above as is usually 

the case for Dirac operators. Hence, the classical variational principle cannot be applied. 
Recently, the Rayleigh functional was used to establish variational principles also for 
Dirac operators in flat space time. The idea is to decompose the given Hilbert space 
on which the Dirac operator is defined into the direct sum of two Hilbert spaces and 
then, in (1), to take the minimum over certain subspaces L in the first Hilbert space, 
and the maximum over certain tp whose first component lies in L. 
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Griesemer and Siedentop [GS99] have proved a variational principle for the eigenvalues 
of a block operator matrix T on a Hilbert space H = H\ ® TL2 in a gap of its essential 
spectrum where the authors did not assume that the operator is semibounded. They 
assumed that the spectra of the operators on the diagonal do not overlap, so that, 
roughly speaking, the decomposition of the Hilbert space Ti into spectral subspaces of 
T is close to the given decomposition. In [GLS99], Griesemer, Lewis and Siedentop 
gave a variational principle for the Dirac operator with Coulomb potential. In the 
case of the block operator matrix in section 4, however, the spectra of the diagonal 
entries do overlap, hence this principle cannot be applied. Under assumptions different 
from those in [GS99], Dolbeault, Esteban and Sere proved a variational principle for the 
eigenvalues of operators in gaps [DESOOa] , in particular they considered Dirac operators 
with Coulomb potential [DESOOb]. Their techniques differ from those used in [GS99] 
and [GLS99], but they also make use of the Rayleigh functional. 

Various types of block operator matrices and their spectral properties have been in- 
vestigated recently. A survey of some recent results can be found in [TreOO]. In this 
paper we are interested in so-called off-diagonally dominant selfadjoint block operator 
matrices 



on a Hilbert space TL = TL\ © H.2 where the linear operators TijiTLj — > TLi) are closed 
and T11 and T22 are symmetric and bounded with respect to Tj* 2 and T12 respectively. 
Further we assume that the operator T\\ is bounded from below by c\ and that T22 is 
bounded from above by 02- However, we do not assume that the spectra of T\\ and T22 
are disjoint. The crucial step to obtain a variational principle for the eigenvalues of T 
is to associate with it an operator valued function, the so-called Schur complement 



The Schur complement plays an important role in the Schur-Frobenius factorisation 
of block operator matrices, see, e.g., [Nag89], [ALMS94], [ALMS96]. It turns out that 
a p (T)n(c2, A e ) = cr p (5i)n(c2, A e ) where A e = inf a ess (Si). To the Schur complement S\ 
we can apply a variational principle proved by by Binding, Eschwe and Langer [BELOO] 
and by Eschwe and Langer [EL04]. We follow an approach that has already been used 
by Langer, Langer and Tretter in [LLT02] for the so-called diagonally dominant case 
with bounded off-diagonal elements T12 and T* 2 . The off-diagonally dominant case has 
been studied also by Kraus, Langer and Tretter [KLT04] under the assumption the 
diagonal entries T\\ and T22 of the block operator matrix are bounded. 
The main result of this paper is a variational characterisation of the eigenvalues of off- 
diagonally dominant block operator matrices T where one of the diagonal entries, e.g. 
T11, of the block operator matrix may be unbounded (Theorem 3.6). Under additional 
assumptions on the operator T12 we derive in Theorem 3.10 the following upper and 





(2) 



Si (A) = Tn - A - Ti 2 (T 22 - A)- 1 T 2 i 



A G p{T 22 ). 



2 



lower bounds for the eigenvalues Ai > A2 > . . . of T which are greater than ci'- 

An < ^Y^^n+no + \] V n +n + \ («21 yfon+no + ||?22|| + Oi) 2 + \ (a + C 2 ) , 
An > V^+no + |( C 1 _ 11^22 1|) 

where the constants 021, a, c\ and C2 are such that 

(y,Tny) < ci ||y|| 2 , y G W 2 , (x, r 22 a;) < c 2 ||x|| 2 , x G Hi, 

||Tnx|| < a||x|| + a2i||?T 2 x||, xeV(T* 2 ), 

and < V\ < 1/2 < . . . are the eigenvalues of T 12 Ti 2 - The index shift no is given by 
no = rmndim£(_ OOj0 )5'l(^) where £(_ aO)0 )5i(A) denotes the spectral subspace of Si (A) 

A>C2 

corresponding to the interval (— 00, 0). In the case T\\ = and T22 = the formulae 
above give the exact values of eigenvalues, hence Theorem 3.12 can be regarded as 
a perturbational result for block operator matrices with a certain type of unbounded 
perturbation. 

The variational principle established in this paper is applied to the angular part of the 
Dirac operator in curved spacetime to derive upper and lower bounds for its eigenvalues 
in Theorem 4.4. These bounds are given explicitly in terms of the physical quantities 
involved. Suffern, Fackerell and Cosgrove [SFC83] derived numerical approximations 
for the eigenvalues by applying a power series ansatz in two of the physical variables 
involved. They obtained numerical approximations for the eigenvalues, but they did 
not give asymptotics of the eigenvalues as are obtained by our analytical approach. 
The paper is organised as follows. In section 2 we consider off-diagonally dominant 
block operator matrices and define the Schur complements associated with them. Sec- 
tion 3 contains the variational principle for block operator matrices which is applied 
to the angular part of the Dirac operator in curved spacetime in section 4. Finally, in 
section 4.3 the analytic bounds are compared to numerical approximations calculated 
by Suffern et al. [SFC83]. 



2 Schur complements and sesquilinear forms related to a 
block operator matrix 

Let Ti\ and H2 be Hilbert spaces with norm and scalar product denoted by || • \\j and 
( " j " )j> J = 1) 2. Consider the Hilbert space H := Ti\ © TL2 equipped with the norm 
|| • || and the scalar product (tp, ip) = (pi, <pi)i + (tp 2 , ^2)2 for tp = (<pi, ^ 2 )*, ip = 
(tpi, t/> 2 ) £ H- On TL we consider block operator matrices 

T = (l 11 1 } 2 ) , V(T) = V x @V 2 C H x ®n 2 

\±2l J 22/ 

with linear operators Tij("Hj — ► Tii), i,j = 1,2. 

In this paper we will always assume that the following conditions concerning the entries 
of T hold: 
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T12 is a closed densely defined operator from H2 to Hi and T-j* 2 = T21; 

^(^12) — ^(^11) an< i ^11 is symmetric in Hi and semibounded from below, 
i.e., there is a constant q £ R such that 

(x,T u x) > ci||x|| 2 , x G £>(Tn); 

T>(T\2) C V(T22) and T22 is symmetric in W2 and semibounded from above, 
i.e., there is a constant C2 G 1 such that 

(x, T 22 x) < c 2 |M| 2 , xeP(T 22 ); 

furthermore, T22 is closed and (02,00) C p(T22). 

We always assume that the block operator matrix T is given by 

(Tl) T=(^ ^2)' V(T)=V(T? 2 )(BV(T 1 2). 

Remark 2.1. (i) Since T\\ is closable by assumption, the condition concerning its 
domain implies (see [Kat80, chap. IV, remark 1.5]) that T\\ is T2i-bounded, i.e., 
that there are positive numbers a and «2i such that 

||Tiix|| < a. \\x\\ + a.21 ||?2ia;||, x G X>(T 2 i). 

(ii) Condition (Dl) implies that T22 is even selfadjoint because the defect index of the 
closed operator T22 is constant on the connected set C \ W(T22) where W(T22) = 
{(x, T22X) : x G V(T22), \\x\\ = 1} is the numerical range of T22. Now, p(T22) D 
C\W (T22) being nonempty implies that T22 has zero defect, hence it is essentially 
selfadjoint. Since T22 is already closed, its selfadjointness is proved. 

Observe that the above conditions do not imply that T is closed. 

Next we associate an operator valued function, the so-called Schur complement, to the 
block operator matrix T . In Corollary 2.8 we show that the spectrum of the Schur 
complement and the spectrum of T are related. For A G pi^Ti'i) define 

V(S [ r n] (X)) := # G V{T* l2 ) : (T 22 - A)" 1 !^ G V(T 12 )}, 

(3) 

5 [mm 1(A) . = Tu _ X _ Tl2 ( T22 _ A) -1 T * 2 . 

The family S"]" 1 (A), A G ^(^22) is the so-called minimal Schur complement of T . In 
the following we are interested in real A. 

Proposition 2.2. Assume that conditions (Bl), (Al), (Dl) and (Tl) hold. Then for 
all A > C2 the operator gj™ 111 ' (A) is bounded from below. If in addition the condition 

(D2) T22 is bounded 
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holds, then (A) is also symmetric, and therefore densely defined and closable. 

Proof. Because of the inclusion D(T^ 2 ) C D(Tn), the Schur complement is well defined. 
To show that 5'{ mml (A) is semibounded, we use that (A — T 22 ) _1 is a positive operator 
for A > c 2 ; hence, for all x 6 T>(Sf am] (X)) we have 

(x, S [ r n] (X)x) = (x, (T n - X)x) - (x, T 12 (T 22 - A) -1 T* 2 x) 

= (x, (Tn - \)x) + (T l2 x, (A - T 22 )~ l T l2 x) > (ci - A)||x|| 2 . 

In particular it follows that the scalar product on the left hand side is real, hence 
5[ mm '(A) is formally symmetric. It remains to be shown that T>(Sf am \X)) is dense in 
Hi. By assumption, the operator (A — T 2 2) _1 is selfadjoint, bounded and positive for 
fixed A > c 2 . Hence there exists a positive square root (X — T 22 )~? which is also bounded 
and selfadjoint. Therefore ((A - T 22 )-^Tf 2 )* = T 12 (X - T 22 )~^ holds. Condition (D2) 
implies that the operator (A — T 22 )~zT* 2 is closed, hence by von Neumann's theorem 
(see, for instance, [Kat80, chap. V, theorem 3.24]) the operator ((A — T 22 )~^Ti 2 )*((X — 
r 22 )~5T 1 * 2 ) = -T l2 (T 22 - Xy l T( 2 with domain 

{xGP((A-T 22 )-lT 1 * 2 ) : (X-T 22 y^ 2 xeV(T 12 (X-T 22 yl)} 

= {xeV(Tf 2 ) : (T^-Xy^x £D(T 12 )} = P(S| min] (A)) 

is selfadjoint and its domain is a core of (T 22 — A)~2Tj 2 ; in particular, its domain is 
dense in Hi. □ 

In order to find selfadjoint extensions of the operators iS[ mml (A) we define sesquilinear 
forms associated with them. 

Proposition 2.3. Assume that the conditions (Bl), (Al) ; (Dl), (D2) and (Tl) hold. 
Then for all X G (c 2 , oo) the sesquilinear form 

V( 5l (X)) := V(T? 2 ), 5l (A)[n, v] := (u, (T u - X)v) - (!>, (T 22 - Xy^v) 

in 7i\ is symmetric, semibounded from below and closed. 

Proof. The symmetry and boundedness from below can be shown as in Proposition 2.2. 
Since T 22 is bounded by assumption, it follows that also (T 22 — A) a is bounded and 
therefore the operator (A — T 22 )~2T 1 * 2 is closed. Consider the auxiliary sesquilinear 
form 

V(t 12 (X)) := X>(2? 2 ), t 12 (X)[u,v] := (!>, (A - T^T^v) . 
For every i\ 2 (A)-convergent sequence (x n .)neN Q T^(T* 2 ) we have 

\\(T 22 -X)~2T* 2 (x n -x m )\\ = k 2 (X)[x n - x m ] 0, n, m ->■ oo. 
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Since (T 22 - X)~^TT 2 is closed, it follows that x := lim x n G V((T 22 - A)^7T 2 ) = 

n— >oo 

V(T* 2 ). This shows that ti 2 (A) is closed. 

The operator Tn — A is symmetric and bounded from below, hence it is form-closable, 
i.e., the symmetric form tn(A) defined by tn(A)[<^, ip] = (Tmp, ip) for tp, ip G T>(Tn) 
is closable; let in (A) denote its closure. Then it follows that the form Si (A) = tn(A) + 
t 12 (X) with domain P(tn(A)) D V(i 12 (X)) is also closed. Since 2?(t 12 (A)) C X>(tn(A)) C 
P(tn(A)), the form Si(A) with domain V(t 12 (X)) = V(T* 2 ) is closed. □ 

Remark 2.4. If we assume that instead of (D2) the condition 

(D2') (T 22 -X)- l (r a n(T* 2 nV(T 12 )) C V(T 12 ), A > c 2 . 

holds, then Proposition 2.2 is also valid. However, the sesquilinear form Si defined in 
Proposition 2.3 is closable but not necessarily closed (cf. [Win06]). 

Throughout the rest of this section we assume that condition (D2) holds. Since the forms 
Si(A) are closed, symmetric and semibounded, there are uniquely defined self adjoint 
operators Si(A) such that V(Si(X)) C D(si(A)) and 

5l (X)[tp, V] = (Si(X)<p, V), <p G ©(5i(A)), V G ^(si(A)). 

In addition, if for fixed ip £ P(si(A)) there exists a w G Hi such that Si(A)[^, = 
(w, ip) for all if belonging to a core of Si(A), then ip € V(Si(X)) and Si(X)ip = w. 
Obviously, for all A G p(T 22 ) the operator Si (A) is an extension of S[ mml (A); it is the so- 
called Friedrichs extension. Note that in general this extension is not the only possible 
selfadjoint extension (cf. [Kat80, chap. VI]). We call the family Si (A), A G p(T 22 ), the 
Schur complement of the block operator matrix T. 

The following proposition gives a sufficient condition such that S[ min '(A) is already 
selfadjoint. 

Proposition 2.5. Assume that the conditions (Bl), (Al), (Dl) ; (D2) and (Tl) hold 
and let A G (c 2 ,oo). In addition assume that condition 

(A2) Tn is symmetric and T± 2 (T 22 — X)~ 1 T^ 2 -bounded with relative bound less than 
1, i.e., there are real numbers a > and 1 > a > (which may depend on 
A) such that for all x G V(T 12 (T 22 - Xy l Tf 2 ) 

\\Tnx\\ < a \\x\\ + 5 \\Ti 2 (T 22 - A) - r x * 2 a;||, 



holds. Then Si (A) = S[ mm ' (A); in particular, S[ mml (A) is selfadjoint. 

Proof. Recall that T\ 2 (T 22 — X)~ l T* 2 is selfadjoint, see proof of Proposition 2.2. If T 22 
is bounded and T\\ is bounded with respect to Ti 2 (T 22 — X)~ l T^ 2 with relative bound 
less than 1, then the Kato-Rellich theorem [Kat80, chap. V, theorem 4.3] yields that 
S| min] (A) is selfadjoint. Hence Sj minl (A) = Si (A) follows. □ 

Next we will show the relation between the spectra of Si (A) and T. 
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Proposition 2.6. For given Hilbert spaces TL\ and TL 2 we consider linear operators 
TijiHj -► Hi), i,j = 1,2, with V(T 21 ) C D(Tn) and P(T 12 ) C V(T 22 ). Let T = 
(tm T22 ) ^ e ^ e ^ oc ^ operator matrix with domain T>(T) = T>(T 2 ±) © T*(Ti 2 ) in the 
Hilbert space H\ © H 2 . If T 22 is bijective, then the operator 

S := T n - T l2 T^T 2l , V(S) := {x G P(T 21 ) : T^T^x G P(T 12 )} 

is we// defined and the following holds: 

(i) T is injective <^=^ 5" is injective 
and 

ker(T) = {(_ r -i T2l/ ) : /Gkersj ker(S). 

(ii) If additionally T 2 \ is surjective, then ran(S) © {0} = ran(T) n (Hi © {0}) and 

T is surjective <4=>- S is surjective. 

Proof, (i) First assume that T is not injective. Then there are / G T)(T 2 i), g G D(Ti 2 ) 
such that 

T n f + T l2 g = 0, T 21 f + T 22 g = and (g)^ ' 

From the second equality it follows that T 22 T 2 \f = —g G T>(Ti 2 ). Consequently, / lies 
in V{S) and / 7^ 0. Inserting the expression for g into the first equality gives Sf = 0, 
hence S is not injective. Now assume that S is not injective and fix an element / 7^ 
in its kernel. For g := —T 22 T 2 \f it follows that 

= 5/ = T u f-T 12 T^T 21 f = T n f + T 12 g, 
= g + T^T 21 f = T 22 \T 22 g + T 2l f). 

Since T 22 is injective, the above equations show / (/, g) 1 G ker(T). 

(ii) For every / G D{S), the element g := —T 22 T 2 \f lies in P(Ti 2 ). Consequently, 
(/, gf G 2?(T) and 



T 



A (T n f + T l2 g\ (T n f -T l2 T 22 l T 21 f\ (Sf 
g) \T 21 f + T 22 g) \ J \ 



which implies that ran(S) © {0} C ran(T) n (Hi © {0}). Conversely, let (/, gf G V{T) 
such that T(/, #)* = (x, 0)* for some x G Hi. From T 2 i/ + T 22 g = it follows that 
g = -T 22 l T 21 f G V{T l2 ). Thus we have / G V(S) and 

x = T n f + T X2 g = Tn/ - T 12 T 22 T 2X f = Sf, 
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implying ran(T) n {TL\ © {0}) C ran(5) © {0}. In particular, the surjectivity of T 
implies that of S. Finally, assume that S is surjective and fix (x, y) 1 G TL\ (&H.2- Since 
ran(T2i) = TL2 by assumption, there is an /' G V{T2\) C D(Tn) such that T^i/' = y. 
Therefore, (/', 0)* lies in the domain of T and we have T(f, 0)* = (Tuf, y) . Since 
we have already shown that ran(5) © {0} = ran(T) n (Hi © {0}), the surjectivity of 
S implies Hi © {0} = ran(T) n (H\ © {0}) C ran(T), hence we finally have (re, yf = 
T(f, 0)' + (re — Tn/', 0)' G ranT because both terms on the right hand side lie in 
ran(T). □ 

The spectrum and resolvent set of an operator valued function are defined as follows. 

Definition 2.7. Let S = (<S(C))c be a family of closed operators, where £ varies in 
some set U C C. Then the spectrum, point spectrum and resolvent set of S are defined 
as 



Analogous definitions apply to the other parts of the spectrum of S, e.g., the essential 
spectrum. 

Recall that for a linear operator S the essential spectrum and discrete spectrum are 
defined by 

a ess (S) := {A G C : dim(ker(S' — A)) = 00 or codim(ran(S' — A)) = 00}, 
o-d(S) := {A G C : A is an isolated eigenvalue of S with finite multiplicity}. 

For a selfadjoint operator S we have crd(S) = \ a ess (S). 

Corollary 2.8. In addition to the assumptions of Proposition 2.2 suppose that the 
operator T is selfadjoint and that T* 2 is surjective. Furthermore, assume that the 
operator function 5| minl defined in (3) is a selfadjoint holomorphic operator family. 
Then we have 



Proof. Proposition 2.6 applied to T — A shows that A G <r(7~) n (02,00) if and only if 
A G cr(S , [ mml ). Moreover, it follows from Proposition 2.6 that A G o~ p (S^ m ^) if and only 
if A G o~ p (T) fl(c2, 00) with dim(ker(T — A)) = 00 if and only if dim(ker(S'[ mml (A))) = 00. 
Hence, (4) is proved. By assumption, for each A G (02,00), the operator ^"""'(A) is 
selfadjoint. Therefore, in order to show (5), it suffices to show crd(T) n (02,00) = 
CT^Sj™ 1111 ) n (02,00). Let A G o~d(T) n (02,00). Then we have dimran(5| mml (A))- L = 
dimker(S'[ mml (A)) = dimker(T — A) < 00. Further, the range of T — A is closed because 



a(S) 
o-p(S) 
P(S) 



{C G U : G (7(5(0)}, 
{C G U : G a p (S(C))} 
{C G U : G p(S(C))}. 



<7 P (T) n p(T 22 ) 
a ess (T) n (c 2 ,oo) 



a ess (5[ mm] ) n (c 2 , 00). 



(4) 
(5) 
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A G ad{T). So Proposition 2.6 shows that ran(S[ minl (A)) = ran(T) n (H\ © {0}) is also 
closed. Hence it follows G <j(i(5[ mml (A)) and consequently A G o"d(5j mml ). 
Let A G <Trf(5[ mml ). Then A G o- p (T) with dimker(T — A) = dimker(5i(A)) < oo and we 
have to show that A is no accumulation point of a(T). Since G crd(S[ minl (A)) and 5[ min] 
is holomorphic, there are 6 > 0, e > and holomorphic functions /Xj : (A — 5, A + 5) — > R 
with //j(A) = for j = 1, ... , dimker(5i(A)), such that for all A G (A — 8, A + 5) we 
have that fj, G cr(S'[ mml (A)) n (— e, e) if and only if // is an eigenvalue of 5[ min '(A) with 
finite multiplicity and /z = /Xj(A) for some j (see [Kat80, chap. IV, §3 and chap. VII]). 
Furthermore, for j = 1, . . . , dimker(S[ minl (A)) we have 

^(A) = ^(x j ,S [ r\X)x J ) = -Wxjf-MTv-Xr^xjf < 

for normalised eigenvectors Xj of S^" 1111 '^) with eigenvalue 0, hence the functions fij are 
not constant in a neighbourhood of A. Consequently, there exists a nonempty interval 
(A — 8, X + 8) such that G p{S l ™ hl] (X)) for all A G (A - 8, X + 8)\ {A}. Consequently, 
<t(T) n (A — 5, A + 5) = {A} which completes the proof. □ 

Remark 2.9. In this section we have considered only the case A G p(T 22 ). For A G 
|o(Tn) the Schur complement 

V{Sf inl (A)) := G V(T l2 ) : (T n - A)" 1 ^ G 2?(7^)}, 
S 2 min] (A) := T 22 - A - r* 2 (r n - A)- 1 ^ 

acting on the Hilbert space H 2 can be used to obtain statements analogous to those 
given for the Schur complements Si (A). 

3 Variational principle for the block operator matrix T 

(T T \ 
rp 1 * Jfa J • 

We assume that the conditions (Bl), (Al), (A2), (Dl), (D2) and (Tl) hold. Let Si (A), 
A G p(T2<z), be the Schur complement of T (cf. Proposition 2.2 and Proposition 2.5). 
The next propositions summarise the properties of the Schur complements and its 
associated forms. 

Proposition 3.1. Consider the selfadjoint block operator matrix 
T =(SS)' V(T)=V{T* 12 )®V{T l2 ) 

on the Hilbert space Hi ®7i 2 . Assume that the conditions (Bl), (Al), (A2), (Dl) and 
(D2) hold. 

(i) For every A G (c 2 , oo), the form 

P( Sl (A)) := V(T? 2 ), 
Si(X)[u, v] := (u, (T n - X)v) - ((!>, (T 22 - A)' 1 ^)) 
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is closed and its domain is independent of X. The operator Si (A) associated with 
the form Si (A) is a well defined selfadjoint operator and S[ mml (A) = Si(A), A G 
(c 2 ,oo). 

Define the operator valued function 

St : (c 2 ,oo) — ^^(Wi), A^5i(A) (6) 
and, for fixed x G 2?(si), i/ie function 

af : (c 2 ,oo) — R, af(A)=*i(A)[x] (7) 
where ^(Tii) denotes the set of all closed operators on 7i\. 

(ii) T/ie operator valued function S\ : (02,00) — > ^(TCi) of (6) is continuous in the 
norm resolvent topology, and for every x G (si) i/ie function af defined in (7) 
is continuous. 

(iii) For every x G D(si)\{0} the function af is decreasing and unbounded from below. 

(iv) If additionally the condition 

(B2) T* 2 is surjective 

is satisfied, then it follows that a ess (Si) = a ess (T) n (c 2 , oo) and <7 p (Si) = a p (T)n 
(c 2 ,oo). 

Proof, (i) The assertions concerning Si (A) are shown in Proposition 2.3 while the iden- 
tity S'[ mml (A) = Si (A) was proved in Proposition 2.5. In particular, the mapping S\ is 
well defined. 

(ii) From (i) it follows that the family of sesquilinear forms (si(X))xe(c 2 ,oo) is of type (a) 
according to the classification in [Kat80]. Hence S\ is a holomorphic family of type (B), 
which implies the holomorphy of Si in the norm resolvent topology. Obviously, for 
every x G ^(si) the function af is even smooth on (c 2 , oo). 

(iii) For every x G T>(si), i / 0, the function af is monotonously decreasing since 

dV f(A) = d~\ Sl(A)N = -ll^ll'-IK^-A)- 1 ^!! 2 < -||x|| 2 < o. (8) 

(iv) This has been shown in Corollary 2.8. □ 

Proposition 3.2. Suppose that in addition to the assumptions of Proposition 3.1 there 
is a constant b > such that for all x G V(T* 2 ) the estimate 

\\n 2 x\\ > b\\x\\ (9) 
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holds. For A G (02, oo) let d(X) be a nonnegative lower bound for (A — T22) 1 , i.e., 

(x, (X-T 22 )- 1 x) > d(\)\\xf > 0, x G W 2 , A G (02,00). (10) 
If there is a 5 > w;zi/i 

5 < d(\)b 2 +ci - A (11) 
/or a// X in a sufficiently small right neighbourhood (c2, C2 + e) 0/C2, i/ien 

(v) i/ie spectral subspace £(_ 00j o)5'i(A) is trivial for all X G (c 2 , C2 + e); 

(vi) ci ess (5i) n (c 2 , c 2 +e) = 0. 

//we a/Zow 5 = in equation (11), i/ien we can show (v) onZy. 

Proof. For A G (02, C2 + e), assumptions (10) and (11) imply for all x G P(5i(A)) \ {0} 
that 

(x,Si(A)x) = Si(A)[x] = (x,T u x)- X\\x\\ 2 + (T* 2 x, (X-T 22 y 1 T* 2 x) 

> (ci - A)||x|| 2 + o!(A)fr 2 ||x|| 2 > <5||x|| 2 . (12) 

(v) If 6 > 0, then for all A G (02, C2 + e) the numerical range of the selfadjoint oper- 
ator Si (A), the closure of which equals the closure of the numerical range of Si (A), is 
contained in the right half plane {z G C : Re(z) > 0}, implying p(5i(A)) 5 ( — 00, 0). 

(vi) If we assume the strict inequality 5 > 0, then the calculation above shows that 
(-oo,5) C p(Si(A)) for A G (c 2 ,c 2 + e), hence (c 2 ,c 2 + e) n<r(5i) =0. □ 

Proposition 3.1 (iii) shows that for every x G ^(si) \ {0} the function o~f has at most 
one zero and is unbounded from below. If in addition (11) holds with some 5 > 0, then 
af is positive for A in a sufficiently small right neighbourhood of c 2 , see (12). Thus the 
continuity of af implies that it has exactly one zero. We denote this zero by p(x), i.e., 

of(A) = ^ X = p(x). (13) 

If relation (11) does not hold, then the function af does not need to have a zero. In this 
case we define p(x) := —00, so that obviously either p(x) = —00 or p(x) > c 2 . Further, 
p(x) does not depend on the norm of x, i.e., for all £ G C \ {0} we have p{x) = p(£x). 

Now fix a linear manifold T> C Tii, independent of A, such that 

P(Si(A)) C V C P(si(A)), AG(c 2 ,co). 

Such a manifold T> exists; for example, we can choose T> = D(T* 2 ). 
For n G N we define the numbers 

a n := min max p(x), (14) 
ice zgLx 
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where L x := L \ {0}. Theorem 3.3 shows that these numbers are indeed well defined. 
Here and in the following, a sequence Ai < A2 < • • • < Atv with N = 00 has to be 
understood as the infinite sequence Ai < A2 < . . . . 

For an interval A C ~R and a selfadjoint operator S we denote its spectral subspace 
corresponding to A by C&(S). By A e we denote the lower bound of the essential 
spectrum of Si, i.e., 



A f 



inf a ess (Si) if a ess (Si) ^ 
00 if a ess (Si) = 



If (C2, A e ) is not empty, then the eigenvalues of Si in this interval are characterised by 
the following minimax principle. 

T* 2 T22 ) with domain T>(T) = 
V(Tl 2 )®V(Ti2) be selfadjoint in the Hilbert space Tti®H.2- Suppose that the conditions 
(Bl), (Al), (A2) ; (Dl) and (D2) are satisfied and that T* 2 is surjective. Further, 
assume that the set (02, A e ) is nonempty and that there is a Ao £ (02, A e ) such that 
dimjC^oo )5'i(Ao) < 00. Then the index shift 

n := mindim/;^ )Si(A) (15) 

A>C2 

is finite and o~(T) n (c2,A e ) consists of a (possibly infinite) sequence of eigenvalues 
Ai < A2 < • • • < Atv where N G No U {00} • If the eigenvalues are counted according to 
their multiplicity, then 

\ n = fi n+no , l<n<N, (16) 
and N S No U {00} is given by 

N = n(A e ) - n Q 
where n(X e ) is the dimension of maximal subspaces of the set 

{x £ V : 3 A > c 2 with Si (A) [re] < 0} U {0}. 



If N = 00, then lim A n = A e . If N < 00 and a ess (Si) = 0, then /i n = 00 for 

n— >oo 

n > no + N. If N < 00, A e < 00, then [i n = A e for n > no + N. 
If there exists a 5 as in Proposition 3.2, then hq = 0. 

Proof. Proposition 3.1 shows that all assumptions of Theorem [EL04, theorem 2.1] are 
satisfied for the Schur complement Si (A), A G (02,00). Hence, the numbers fJ> n +n exist 
and are equal to the eigenvalues of the operator family Si. By Corollary 2.8, we have 
o-p(Si) = o-p(T) n (c2,oo) and a ess (Si) = o ess (T) n (02,00) so that all the assertions 
follow from theorem [EL04, theorem 2.1]. 

If even the assumptions of Proposition 3.2 are valid, then it follows automatically that 
(c2j A e ) 7^ and that dim J C(_ OO)0 )<S'i(A) = for A in a sufficiently small right neighbour- 
hood of C2, hence the index offset no appearing in formula (16) vanishes. □ 
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The numbers p(x) are rather hard to estimate. However, there is a representation of 
p{x) as the supremum of a functional A+ (y) where y varies in some subspace of Ti.2, 
see (17) and Lemma 3.5. The functional A+ is connected with the so-called quadratic 
numerical range of block operator matrices, see, for example, [LT98] and [LMMT01]. It 
was used in [LLT02] to obtain a variational principle for block operator matrices with 
bounded off-diagonal entries. 

(rp rp \ 
T2 T22 ) be a closed block operator matrix on the Hilbert 

space Ti = Hi © H2 with domain T>(T) = T>{T2\) © T)(T\2). Assume that the operators 
T\2 and T21 are closed and that Tu is T2\-bounded and that T22 is T\2-boundcd. For 
(x, y) 1 G T>(T), 1 / 0, 1; / 0, consider the matrices 

' (x,T llX ) (xjTnyy 



N \\y\\ 

y, T221 

IMI llvll llyll 2 



T xv := G M 2 (C) 

,V 1 (y,T 2 ix) (y,T 22 y) 1 



with eigenvalues 



, x\ 1 / (x, T xl x) (y, T 22 y) 
±{ y) : " ^HRF" HP* - 



± I i (x,T 1} x) (y,T 22 y) \ 2 ! 4(x,T 12 y)(y,T 21 x) \ ^ 



and define the sets 



A±(T) := U±(*) : x G P(T 21 ) X , y G ©(T, 



12, 



L X 



where we use again the notation L x := L\{0} for linear spaces L. Note that in general 
the radicand may be negative or complex. The quadratic numerical range W 2 {T) ofT 
is defined as the set of all complex numbers A that are eigenvalues of some is. 

W\T) := |J a p (T x>y ) = A + (T)UA_(T). 

xSX>C T 2l) X 
j/6B(T 12 )X 

The number A± ( y ) does not depend on the norm of the vectors x and y. Therefore 
it suffices to restrict the definition of A± (y) to elements (x, y)' G £>(T) with ||x|| = 

l|y|| = i- 

In the following we characterise p(x), defined in (13), in terms of A± (y). Recall that 
p(x) is the unique zero of the function A h- > crf(X) if it exists and p{x) = — oo otherwise. 

Lemma 3.5. Assume that the conditions of Proposition 3.1 hold. Then for all x G 
P(T* 2 ) \ {0} with p{x) ^ — oo we have 



p{x) = su P {a + Q : yGP(T 12 )\{0}}. 



(18) 
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If in addition x G V(Si(p(x))), then the supremum is attained, thus we have 

p(x) = max{A + Q : y G V(T 12 ) \ {0}}. (19) 

Proof. Fix x G T>(T* 2 ) \ {0}. Since T 2 i = T* 2 and the operators T\\ and T 22 are 
symmetric, (17) shows that \ + {y) is real for all y G T)(Ti 2 ) \ {0}. Note that for 
arbitrary A > c 2 , x G P(T* 2 ) \ {0}, y G P(T i2 ) \ {0} we have 



|x|| 2 ||y|| 2 det(7i i3/ - A) 



= ((x, Tux) - A||x|| 2 )((y, T 22 y) - A||y|| 2 ) - (x, T 12 y)(y, T£ 2 x) 

To prove the assertion we first show that p(x) > A+ (y) for all y G V(Ti 2 ) \ {0}. So 
fix y G V{T\ 2 ) \ {0} and, for simplicity of notation, set A + := X + ('y). If A + < c 2 , then 
nothing has to be shown since c 2 < p(x) by assumption. Now assume A + > c 2 . Since 
p{x) is the unique zero of the monotonously decreasing function af, it suffices to show 
o~f (A+) = Si(A+)[x] > 0. By definition, A+ is an eigenvalue of the complex 2 x 2-matrix 
T x , y , thus, by (20) and the definition of Si(A) in Proposition 2.3: 

= ||x|| 2 ||y|| 2 det(T x>y - A+) 

= (y, (T 22 -X + )y) 5l (\ + )[x] (21) 
+ (y, (T 22 - A+)y) (T\* 2 x, (T 22 - A + )^ 1 T 1 * 2 x) - |(y, T^x)) 2 . (22) 

For A > c 2 the operator (A — T 22 ) is strictly positive and the same holds for the induced 
sesquilinear form (u, v) 1— > (ti, (A — T 22 )u) for u, v G T>(T 22 ). For this form we have the 
following generalised Cauchy-Schwarz inequality 



(«, (A - T 22 )v) | 2 = I ((A - T 22 )^u, (A - T 22 )^v) 



2 



< ||(A-T 22 )5n|| 2 ||(A-T 22 )5t;|| 2 
= (u, (A - T 22 )u) (u, (A - T 22 )i>) 

for all u, v G V(T 22 ). Since y G D(Ti 2 ) C "D(T 22 ), we can use this inequality to estimate 
the two terms in (22): 

(y, (T 22 - \+)y) (T^x, (T 22 - A + )~ 1 T 1 * 2 x) - |(y, T^x)) 2 

= (y, (T 22 - A+)y) (T^x, (T 22 - A+)^ 1 T 1 * 2 x) 

- I (y, (A+ - T 22 )(A+ - T 22 )- 1 T 1 * 2 x)| 2 

> (y, (A+ - T 22 )y) (T^x, (A+ - T 22 )^ 1 T 1 * 2 x) 

- (y, (A+ - T 22 )y) (T^x, (A+ - T 22 ) _1 T* 2 x) 
= 0. 

Because the factor (y, (T 22 — A + )y) in the term (21) is negative, it follows that the 
second factor, si(A + )[x] = af (A + ), must be nonnegative, and thus we have proved the 
inequality p{x) > sup{A + ( y ) : y G V(Ti 2 ) \ {0}}. 
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If x G V(Si(p(x))), then we can choose an element y such that p(x) = X + (y). To this 
end, define y := (T 22 —p(x)) T± 2 x. This vector is well defined and it lies in the domain 
of T\2 since by assumption x E V(Si(p(x))). If we use 

(y, (T 22 -p(x))y)(T^x, (T 22 - pix^T^x) - \{T? 2 x, y)\ 2 = (23) 
and 5i(p(x))[x] = 0, it follows with the help of (20) 
||x|| 2 ||y|| 2 det(7^ > 2 / - p{x)) 

= {y, (T 22 - p(x))y) Si(p(x))[x] 



+ (y, (T 22 -p(x))y){T* l2 x, (T 22 - p{x))~ l T* l2 x) - \(y, T{ 2 
0. 



x)\ 2 



This implies that p(x) is an eigenvalue of T x ,y Together with A_ ( % ) < A + ( % ) < p(x) 
it follows that p{x) = \ + {y) which proves (18) and (19) in the case x E V{S\(p{x))). 
It remains to show (18) in the case x ^ T>(Si(p(x))), i.e., for elements x E D(T^ 2 ) such 
that (T 22 — p(x))~ 1 T^ 2 x ^ V{Ti 2 ). If x E V(T^ 2 ) \ V(Si(p(x))), there exists a sequence 
(x n ) nS M C V(Si(p(x))) such that 

x n — ► x and (T 22 - A)^2T 1 * 2 x n — > (T 22 - A^T^x, n — ► 00, 

since in the proof of Proposition 2.2 we saw that T>(Si(p(x))) is a core of (T 22 — 
p(x))-5Tf 2 . Set y n := (T 22 - p(x))" 1 T 1 * 2 x n , n G N. Because both (T 22 - p(x))^ 1 
and T 22 — p(x) are bounded, the limites y := lim y n and lim T? 2 x n exist and are not 

n— >oo 

zero; otherwise it would follow that x G D(Si(p(x))) in contradiction to the assump- 
tion on x. Moreover, since T\\ is relatively bounded with respect to T* 2 , also the limit 
lim T\\x n exists. Therefore, all terms in 

n— >oo 

Si(p{x))[x n ] 

= Si(p(x))[x] +si(p(x))[x n - x,x n ] +si(p(x))[x n ,x n - x] - 5i(p(x))[x n - x] 

= 2 (x n - x, (Tn -p(x))x n ) - (x n - x, (Tn -p(x)){x n - x)) 
+ 2 (Tl 2 (x n - x), (T 22 - pix^T^Xn) 
- { T 12( x n - x), (T 22 - p(x))~ l Tl 2 (x n - x)) 

converge to zero for n — ► 00. As in (23), we obtain 

(y n , (T 22 - p{x))y n ) (T* 2 x n , (T 22 - p{x))- l T{ 2 x n ) - |(T* 2 x„, y n )\ 2 = 0, tieN, 

which implies 

IMI 2 \\y n \\ 2 det(T Xn: y n -p{x)) = (y n , (T 22 - p(x))y n ) Si{p{x))[x n ] n ^°° > 0. 
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Since neither x n nor y n tend to zero, it follows that 

(p(x) - A_ (j))(p(x) - A + = det(T Wn -p(x)) -=^=_> 0. (24) 

Each entry of 

T - T 
2 x n ,yn 2 x,y 

(x n , Tn{x n - x)) + (x n - x, Tux) (x n , T 12 (y n - y)) + (x n - x, T 12 y)\ 

(y n , T^ 2 (x n - x)) + (y n - y, T£ 2 x) (y n , T 22 (y n - y)) + (y n - V, T 22 y) J 

(x n , Tu(x n - x)) + (x n - x, Tux) (Tf 2 x n , y n - y) + (Tf 2 (x n - x), y)\ 
(yn, Tf 2 (x n - x)) + (y n - y, T% 2 x) (y n , T 22 (y n - y)) + (y n - y, T 22 y) J 

converges to zero for n — > oo, hence we have T Xn ,y n — ► T xy in norm. Thus the eigenvalues 
A-t (y^} °f ^in,yn converge to the eigenvalues X± ^\ of T x>y , in particular it follows 
that 

p(x) = lim A_ 



n— >oo 



Since X± is continuous in both its independent variables, and since y n £ T>(T^ 2 , it 
follows that 

p(x)= Km A+(*) <sup{A + (*) :^^)\{0}}. (25) 

□ 

The following theorem is the main theorem of this paper. It provides a variational 
characterisation of the eigenvalues of T in a right half plane of C in terms of its entries 

Theorem 3.6. Suppose that the assumptions of Theorem 3.3 hold, that is, suppose that 
conditions (Ti), (Bl), (Al), (A2), (Dl) and (D2) are fulfilled, that T* 2 is surjective, 
that (c 2 , A e ) 7^ and that there is a Xq G (c 2 ,X e ) such that dim£(_ 00) o)<S'i(Ao) < oo. 
Then the eigenvalues ofT in (c 2 ,X e ) are given by 

X n = min max sup A + ( X ) , 1 < n < A, (26) 



ICC(T« 2 ) i£Lx y ev(T 12 ) x \y 

dim L — n-\-TiQ 

where we have adopted the notation of Theorem 3.3. If the domain of Si (A) does not 
depend on X, i.e., if 

V(S l (X)) = V(S 1 ), AG (c 2 ,oo) 

then we have 

X n = min max max A+ ( X | , 1 < n < A. (27) 

LCOfSi) x gLX y<zV{T 12 ) x \V ) 

dim L — n+720 
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Proof. From Theorem 3.3 it follows that all eigenvalues of T greater than C2 are given 
by 



A n = min max p(x), 1 < n < N, 

ICC 

dim L— ri + Ti,Q 

where T> is any linear manifold with 2?(Si(A)) CDC T>{s\). From Proposition 3.1 we 
know that the forms Si(A), A G (c2,oo), are closed and that £>(si(A)) = P(T 1 * 2 ). Fix 
n > and a subspace L C 2?(si) with dimL = n + n,Q. Then there exists an x G L with 
p(x) ^ — oo. Lemma 3.5 yields 

max p{x) = max p{x) 

= max sup A + ( X J = max sup A + ( 

PfL* y£T>(T 12 )x \V J zSLx (T )x 

If we have even L C P(5i), then the supremum can be replaced by the maximum. □ 

In general, it is not easy to determine the index shift hq. Sufficient conditions for the 
finiteness no, which are met the operator in the application in Section 4, are given in 
the following proposition. 

Proposition 3.7. (i) If in addition to the assumptions in Theorem 3.6 the operator 
T^ 2 is compactly invertible and the operator T\\ is bounded, then the index shift 
no is finite. 

(ii) If the assumptions of Proposition 3.2 are satisfied, then no = 0. 
Proof. Let A G (c2,oo). 

(i) The assumptions imply that T\2i^22 — ^) -1 7\*2 ^ s compactly invertible. Since T\\ is 
bounded, also Si (A) is compactly invertible, hence its spectrum consists of a sequence of 
eigenvalues with finite multiplicity which has no accumulation point. Since the operator 
S'i(A) is bounded from below, it follows that dimC^ oo ^Si(X) < oo, in particular, 

no = min dim C(_ 00 )Si(A) < oo. 

A>c 2 

(ii) Proposition 3.2 implies that £(_ OOj0 )5'i(A) = for A sufficiently close to C2, hence 
n = min A>C2 dim £( oo 0) Si (A) = 0. □ 

In the application in section 4, the spectrum of T±2T* 2 consists of simple discrete eigen- 
values only. For this situation, we specialise Theorem 3.6 further. 

Remark 3.8. Assume that T = ( T ° 2 T q J with domain V{%) = D(T* 2 ) © V(T 12 ) C 

Hi © Ti-2 is closed and that T\2T* 2 and T* 2 T\2 are strictly positive. Then o~ p (To) = 
{A G E : A 2 G a p (T 12 T? 2 )}. 
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Proof. For A £ o" p (7o) \ {0} we have A 2 £ ^(T^T-J^) n ^(T^T^) since for each eigen- 
vector (/, gf of T with eigenvalue A it follows that / £ V(T 12 Tf 2 ), g £ V(T? 2 Ti 2 ), 

f, g + and = (T + A) (T - A) ( { ) = ( ( ^ ^ ) - A 2 ) ( J ) = ( ) . 

On the other hand, if \i 7^ is an eigenvalue of T^ 2 T\ 2 with eigenfunction g, then it is also 
an eigenvalue of T\ 2 T\ 2 with eigenfunction T\ 2 g. For cr = ±1 we define / = o~[i~2 T\ 2 g. 
Then we have that (7o — (Jy/Ji) (£) =0, hence are eigenvalues of Tq. □ 

To estimate the functionals A+ ( y ) , we use the following auxiliary lemma. 

Lemma 3.9. For a\, a 2 , b\, b 2 , 7 £ R wrai/i ai < 61 and 02 < 62 we define the function 



f : [01, 61] x [a 2 , 62] R, /(s, t) = s + t + ^(s-O' + T 2 - 
For fixed t, the function f is monotonously increasing in s and vice versa. In particular, 

/(ai,a 2 ) < f(s,t) < /(6i,6 2 ), (a, t) £ [at, h] x [a 2 , 62]. 
Proof. Partial differentiation of / with respect to s yields 



= i + ^^= > v^^a^ > 0. □ 

Theorem 3.10. Lei T = ^2) with domain v { T i2) © ^(^1*2) QHx®H 2 be a 
self adjoint block operator matrix such that the conditions (Tl), (Bl), (Al), (A2), (Dl) 
and (D2) hold. Then T\\ is bounded with respect to T* 2 ; let a, a 2 \, c 2 and c 2 such that 

\\T n x\\ < a\\x\\ + a 2 i\\T* 2 x\\, xeV(T* 2 ), 

c 2 \\y\\ 2 < (y, T 22 y) < c 2 \\y\\ 2 , y £ H 2 . 

Further, let T^ 2 be bijective and assume that for all A £ (c 2 , 00) the Schur complement 

Z>(Si(A)) = {xeV(Tf 2 ) : (T 22 - A)" 1 ^ £ V{T 12 )}, 
Si(X) = Tii-A-TiaCZfea-A)- 1 ^, 

is selfadjoint and that T>(S\(X)) =: "D(Si) is independent of A. Additionally suppose 
that there exists a Xq £ (c 2 ,oo) such that dim£(_ oo )5i(Ao) < 00. // the spectrum of 
the operator Ti 2 T* 2 satisfies 

a(T 12 T* 2 ) = a p (T 12 T* 2 ) = { Vj : j £ N} with < v x < v 2 < . . . 

where the eigenvalues are counted with their multiplicities, then the block operator matrix 
T has discrete point spectrum X± < X 2 < . . . A^r in (c 2 ,X e ). More precisely, if no is as 

in Theorem 3.3, that is, uq = min dim£(_ 00 )5i(A), then for all 1 < n < N the 

A>c 2 ' 

eigenvalues X n ofT in (c 2 ,X e ) satisfy the estimates 

An < ^-V^n+no + J ^n+n + I («21 y^n+no + |« ~ C 2 |) 2 + \ {p. + C 2 ) , (28) 



An > A/ ^n+no + \{c\ ~ C 2 ) 2 + \{ci + C 2 ). (29) 
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Proof. Since T* 2 is closed, its resolvent is bounded by the closed graph theorem. Hence 
Proposition 3.2 yields that (c2,A e ) 7^ so that all assumptions of Theorem 3.6 are 
satisfied. In particular, the index shift no is finite. To prove inequalities (28) and (29), 
we estimate the right hand side of (27). Note that T>(T 22 ) = W 2 arid that 

(i) (x, Tux) < \(x,Tnx)\ < \\x\\ \\Tnx\\ < \\x\\ (a\\x\\+a2i\\T^ 2 x\\), ieP(Tj 2 ), 

(ii) (x, Tux) > ci||x|| 2 , x G V(Tn), 

(iii) (y, T 22 y) < c 2 ||y|| 2 , y G H 2 , 

(iv) (y,T 22 y) > c^" 1 1 2/ 1 1 2 , 

(v) \(y, T^ 2 x)\ 2 < \\y\\ 2 \\T* l2 x\\ 2 , x G V(T* 2 ), y G H 2 . 

First we prove (28). With the help of inequalities (i), (iii) and (v) and the auxiliary 
Lemma 3.9 we find for x G V(T? 2 ) \ {0}, y G V(T 12 ) \ {0}: 



x 



a + i~ \ = - r , + v „ ,tr y +w v , ' - ^ ;r y + 



+ * y J 2 1 ||x|| 2 ||y|| 2 y \ ||^|| 2 ||y|| 2 / \\ x \\ 2 IMI 2 



1 / (x, Tux) (y, T 22 y) / / (x, Tux) 

(y, T22V)\ 2 , 4|(y, t; 2 x" 2 



o I ' II II 1 z 1 W \ 1 II II z I 1 II 119 II 119 

2 1 ||x|| y \ ||x|| y IfII II2/II 



„ 1 / a 2 i||T 1 * 2 x|| //a 2 i ||T*x|| A 2 4||T*x|| 2 

2 I ||x|| y \ \\x\\ J \\x\\ 

The right hand side is independent of y and monotonously increasing in ||T* 2 x||. For 
given n G N let «Sf n be an n-dimensional subspace of the spectral space Ly vi jl/n ](Ti 2 T^ 2 ) . 
Then for every x G Jz? n we have that HT-^xH 2 = (x, T\ 2 T* 2 x) < v n \\x\\ 2 . 
Observe that T){si) = 2?(Ti 2 )*, thus the minimax principle (27) shows that 



X n = min max sup A+ ( ) < max sup A+ 



dim L— n + nQ n + nQ 



X 



I a 2 i\\T? 2 x\\ /a 2 i\\T*x\\ A 2 4||T*x|| 2 



< max -I Q + c 2 + +1/ n 2 +\a-c 2 \) + 



2 



c cv> x 2 \ x V V x / x 

< i + C 2 + a 2 l7^n+n + y/ (a21\AVfn + |« — C 2 |) 2 + 4 f n+rlQ J 

which proves (28). Note that if o^i^/^no > a — c 2 , then the above estimates is true with 
a — c 2 instead of \a — c 2 \. 

In order to show (29), we choose a particular y G V{T\ 2 ). Since by assumption T^ 2 l 
exists and is bounded by b~ 1 , also T^ 2 exists and is bounded by ft -1 . For every x G 
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T>(T* 2 ) the element y{x) := T 12 x exists and lies in T>(T\2)- Therefore, again by (27) 
and the inequalities (ii) and (iv) we obtain 



A n = min max sup A_ 

LQT>(T* 2 ) l6 LX y£V{T 12 ) x 
dim L — n+TiQ 



> mm max sup - I c\ + c 2 \ / (ci — c 2 J z H r — 2 2 — ) (30) 

lc-d(t* 2 ) x <=Lx yev(T 12 ) x ^ V v IfII Hj/II 



dim Zj — n+n 
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1/ /. _ N9 4 fx, x) 2 \ 

> mm max - ci + c 2 + W (ci - c 2 ) 2 + — I 

ic©(T* 2 ) x&l x 2\ y ||T 12 x|| 2 \\x\\ z J 

dim L — n+TiQ 



1 



mm max - I c\ + c% + y (c\ — c 2 ) 2 + 4 ||T 12 1 x|| 2 ||x" 2 

iCD(T* 2 ) XGL X 2 
dim L — n+rag 

= ^(ci + c 2 ) + \ \{ci - c 2 ) 2 + C min max IIT^ 1 ^!!" 1 \\x\\) . (31) 

2 V ^ LC-D(T*) x£L x " J 

dim L — n+TiQ 

For every n-dimensional subspace L n C 2?(T* 2 ), also L n C T>(T* 2 Ti2) is n-dimen- 
sional. Hence it follows that 

min max 1 1 ^"xa 1 1 1 1 ll^ll = mm max 11^12^11 

LQT>(T* 2 ) X&L X LCT)(T* 2 ) ^T~ 2 L X 

dim L=7i-\-riQ dim L — ti+tiq 

= min max ||Ti 2 £||. 

LCT>(T* 2 T 12 ) ggLX 
dim L=n-\-riQ 

The squares of the nonzero eigenvalues of To = ( ^ T * 2 ^ = ^ °j q ^ 7o ^ ^ o ) are ^ ne 
eigenvalues z^i < vi < . . . of T^T^ (see Remark 3.8). On the other hand, the variational 
principle of Theorem 3.6 applied to Tq shows that 

r- > • \(y 1 TuO\ 
Jv n = \ n = mm max max — - — — 

LCT>(T* 2 T 12 ) 5eL x yeV(T* 2 ) x \\y\\ ||£|| 
dim L — n 

\(T 12 £, T 12 Q\ . \\T 12 g 

< mm max — rnrr, — = mm max — — - — . 

LeX.(T* 2 T 12 ) ?eL x \\T 12 $\\ U\\ LeT>(T h T 12 ) ?eL x ||^|| 

dim L — n dim L — n 

Inserting into (31) yields 



An > \{ci + C 2 ) + J \(c\ - C 2 ) 2 + U n+no . □ 

Remark 3.11. Theorem 3.10 can be regarded as a perturbation result for the eigenval- 
ues of the block operator matrix y j under the unbounded perturbation ( T ^ ^ 

since in the case T\\ = T22 = the spectral shift uq vanishes and the estimates (28) 
and (29) reduce to A n = yjv^,. If the sequence (v n ) n of the eigenvalues of T±2T* 2 is 
unbounded, then A has the same asymptotics as \/u^, i.e., . Xn — ► 1 for n — ► 00. 
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If also the operator T\\ is bounded, then the estimate for \ n from above can be further 
improved. 

Theorem 3.12. In addition to the assumptions in Theorem 3.10, let Tn and T 22 be 
bounded. Then there are real numbers c\ and cf such that 

ci ||x|| 2 < (x, T\\x) < ||x|| 2 , x G 7i\. 

Let no = mindim£/_ 00 ")5i(A). Then the eigenvalues of the block operator matrix T 

A>c 2 v ' ' 

in (c 2 , A e ), enumerated such that c 2 < X\ < A2 < . . . , can be estimated by 

An < \J Vn+n + l{cf - C 2 ) 2 + ^(c+ + C 2 ), 1 < Tl < N, (32) 

An > ^„+„ + |(ci -C2) 2 + |(ci + cj), 1 < n < N. (33) 

where < 1^1 < v 2 < . . . are the eigenvalues ofT 12 T^ 2 , see Theorem 3.10. The index 

shift uq is given by hq = min dimjC^^ )5'i(A). 

A>c 2 v ' ; 

Proof. We only need to show (32). If we use (x, T\\x) < cf \\x\\ 2 , x G TL\ instead of 

inequality (i) in the proof of formula (28), we obtain with the help of the auxiliary 
Lemma 3.9 

A+ ( I ) = \ TllX) + (y ' T22V) + ^ {{X ' TllX) ~ T22y))2 + 4|(?/ ' Tl * 2X)|2 ) 
< ^c+ + c 2 ) + ^/i( C +- C2 ) 2 + ||T 1 * 2 x|| 2 

for all (x, yY £ T>(T) with = ||y|| = 1. Now formula (32) follows by a reasoning 
analogous to that of the proof of (28). □ 



4 Application to the angular part of the Dirac equation 
in the Kerr-Newman background metric 

The Kerr-Newman metric describes the spacetime in the exterior of an electrically 
charged rotating massive black hole. A spin-| particle with mass m and electrical 
charge e outside the black hole obeys the Dirac equation 

(§ + <K)$ = (34) 

where ^ is a four-component wave function describing the particle and 21 and 9\ are 4x4 
differential expressions that contain partial derivatives with respect to all four spacetime 
coordinates. It can be shown that by a suitable ansatz the Dirac equation (34) can 
be decoupled into a system of two ordinary differential equations ([Cha98], [Win06], 
[WY06]): the radial equation that contains only derivatives with respect to the radial 
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coordinate and the angular equation that contains only derivatives with respect to the 
angular coordinate For recent results on the radial equation see [Sch04] and [WY06]. 
The angular equation is given by 



(2l-A)* = on (0,tt). 



with the differential expression 



21 = 



( 



am cos ■!? 



am cos $ 



) 



on (0,tt), 



(35) 



where 



d k + 1 



+ auo sin 



(36) 



d# sin 



The number k G Z describes the motion of the electron in the plane of symmetry. The 
parameter a := J/M G M describes the rotation of the black hole where J is the angular 
momentum and M is the mass of the black hole. 

For the following results on the angular operator we refer to [Win06]. Let 



with the norm \\(f, g) l \\ = yjj/jh + lblb where || • ||2 denotes the usual norm on 
J2? 2 ((0, 7r), d$). It can be shown that on ft the formal expression 21 has the unique 
selfadjoint realisation 



where D is the operator of multiplication by the function (0,7r) — > M, ~& i— ► am cos i? 
and -B is the first order differential operator given by 



B is closed and its inverse can be computed explicitely. It turns out that B is a 
Hilbert-Schmidt operator, hence it is compact. Consequently, also the angular operator 
A in the special case a = is compactly invertible. Since all terms in A involving 
the parameter a are bounded, a perturbation argument shows that A is compactly 
invertible for every value of a. Furthermore, it can be shown that the spectrum of A 
consists of simple eigenvalues only. 



Remark 4.1. It can be shown that if i? i— > (/(??), 3(1?))* is an eigenfunction of A, then it 
follows that there is a 7 G C with (7! = 1 such that (c/(vr-??), f(ir-'&)) t = 7(/(tf), 5(1?))*, 
1? G (0,tt). 

In the special case a = the eigenvalues can be calculated explicitly; one obtains 
a p (A) = {sign(n)(|A; + \\ - \ + n) : n G Z \ {0}} for a = 0. 



ft := J^ 2 ((0,vr),dtf) x J^ 2 ((0,vr),dtf) 





(37) 



V{B) = {/ G ft : / is absolutely continuous, »+/ G ft}, = »+/. 
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The operator BB* is a Sturm-Liouville operator with spectrum consisting of discrete 
simple eigenvalues < v\ < U2 < • • • only. By Sturm's comparison theorem we obtain 
the following two-sided estimates for the eigenvalues v n of T^T^'- 

max{0, (\k + \\ - \ + n) 2 + < v n < (\k + \\ - \ + n) 2 + Q + (38) 

with 

[a 2 J 1 + \ + 2(k + \)ouj if 2au 4. [-1, 1 ], 
Q_ = 2(fc + |)au;-| aw l> 0+ = < (39) 

I 2(/c + ~)aw + \au)\ if 2au G [—1, 1 ]. 

Remark 4.2. In the case a = these estimates give the correct eigenvalues of BB*. 
The corresponding eigenfunctions are hyper geometric functions. 

For the block operator matrix A, the Schur complement for A £ p{D) is given by 



©OSi(A)) = {/£ : (-D — \)~ x B*f G P(B)}, (40) 

5i(A) = -£>- \-B{D- X)-^*. (41) 

Lemma 4.3. T/ie angular operator fulfils conditions {Tl), (Bl), (B2), (Al), (A2), 
(Dl) and (D2) of the preceding section, in particular, we have 

c± := — \am\ < (x, —Dx) < \am\ =: cf , x £?i (Al') 

C2 := — |am| < (x, Dx) < \am\ =: C2, x G TC, (Dl ; ) 

= =\am\, (A2'), (D2') 



and <r(D) = a{—D) = a ess {D) = [— \am\, \am\ ]. For all A G (|ani|, 00), the form 

V( 5l (X))=V(B*), B 1 (X)[f,g] := (/, (-D — X)g) — (B*f, (D - X)~ 1 B*g), (42) 

is symmetric, semibounded from below and closed. Further, the operator Si (A) is the 
self adjoint operator associated with Si(A) ; and its domain is independent of X, more 
precisely, we have 

V{Sx(X)) = V(BB*), AG(|om|,oo). (43) 

Proof. Since D is the bounded operator given by multiplication with the continuous, 
nowhere constant function amcos-d, $ G (0,7r), the assertions concerning the spectrum 
of D and relations (Al'), (Dl') and (A2'), (D2') are clear. Hence conditions (Al) and 
(Dl) are satisfied with C2 = \am\ and c\ = — \am\, and (A2) and (D2) hold because 
D is bounded with \\D\\ = \am\. Since a(D) = [—\am\,\am\], the sesquilinear forms 
Si(A), A G (|am|,oo), are well defined, and, by Proposition 2.3, they are symmetric, 
semibounded from below and closed. Proposition 2.5 implies that for A G (|am|, 00) the 
operator Si (A) is the selfadjoint operator associated with Si(A). 
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To prove (43), fix / € V(BB*) and A £ (|om|,oo). We have to show (D - A)~ 1 J B*/ £ 
T>(B). Since both (am cos 1 !? — A) -1 and B* f are absolutely continuous, we have 



%$+(d — \y 1 B*f($) 

(d A; + — \ / d k + — \ 
— H j- + awsin-i? ) (am cos ■& - A) -1 ( — — H f + au;sini9 )f0&) 
du sinw / \ d» sinw / 

(d k + — \/d + — \ 
— H f + awsintf ) ( — — H f + awsm??)/(i9) 
dv smv / \ di/ sun? / 

+ (^ (ara< „*-A)-')(-± + i±i + „„sin*)/ M 

= (I? - A) _1 5B*/W + (^(amcostf- A)- 1 )^*/^)- 

Observe that the first term on the first line is the formal differential expression as- 
sociated with B. Since, by assumption, / £ T>(BB*) and since both (D — A) -1 and 
g^(am cos $ — A) -1 are bounded operators on H, it follows that (D — X)~ X B* f £ T>(B), 
and consequently / £ T>(S±(X)). 

Conversely, assume / £ T>(Si(X)) for some A £ (|am|, oo). The function am cos $ — A is 
differentiable on (0, tt), hence we have 

f ^ k -\- - \ 

( — H — : — ^ + au> sini? ) (am cos i? — A)(amcos i? — A) -1 



V di? sin ■& 

V an smv J 

/ d /c + - \ 

(am cos i? — A) ( — H ^ + aw sin $ ) (am cos i9 — A) -1 

V di? sin v J 

-fr + ^ + au, sin i?)/(i?) 
di/ sin u / 

+ f~r7T( am cos ^ ~~ ■^) N ) ( am cos $ ~~ f — wT "I + aw sini?) /(i?) 

\dt7 / V di; sun? / 

= (D - X)B{D - \)- l B*f($) + (-(arocosi?- A))(£> - X)~ 1 B*f{v'). 

Since D — X and (am cos i? — A) are bounded operators on it follows that the 
function above is also an element of TC, hence we have B*f £ T>(B) which implies that 
feV(BB*). □ 

4.1 Explicit bounds for the eigenvalues of A 

As mentioned already earlier the spectrum of the angular operator consists only of 
isolated simple eigenvalues without accumulation points in (—00,00). We also know 
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that the eigenvalues depend continuously on the parameter a. To express this depen- 
dence explicitly we frequently write A n (a). Therefore we can enumerate the eigenvalues 
X n (a), n £ Z\ {0}, unambiguously by requiring that X n (a) is the analytic continuation 
of A n (0) = sign(n)(|/c + ^ | — ^ + |n|) in the case a = 0. Since all eigenvalues are simple, 
it follows that A n (a) < A m (a) for n < m. 
For fixed Kerr parameter a we define m± £ Z such that 



• •• < A m __ 2 (a) < A m __i(a) < -\am\ < A m _(a) < ... 

• < A m+ (a) < \am\ < X m++1 (a) < A m++2 (a) < ... 

i.e., cr(A) n [— |am|,|am|] = {A n (a) : m_ < n < m + , n ^ 0} and the number of 
eigenvalues of .A in the interval [— \am\, \am\ ] is given by 



# (o"(-4.) n [— |am|, \am\ ] ) 



m + — m_ if € [m_ , m + ] 

m_|_ — m_ + 1 if ^ [m_ , m + ] 



Observe that m + and m_ depend on the physical parameters a, m, uj and /c. 

Theorem 4.4. Let uq = dim£(_ oo )5'i(Ao) for some Ao £ (|am|, A m++ i). T/ien the 
eigenvalues of the angular operator A to the right of \am\ are given by 

A mj _+n(a) = rnfn max max A+ ( X 1 , n £ N. (44) 

+ LC-D(BB-) x£L x y£V(B)x \yj 

dim L— n+riQ 

Furthermore, the eigenvalues can be estimated by 

^Vn Q +n ~ \am\ < X m++n (a) < y/v no+n + \am\, n £ N, (45) 

where v n + no are the eigenvalues of BB* . Explicit estimates for \ n in terms of the 
physical parameters a, m and oo are 



with 



a!! ] + „ < A m++n (a) < X [ :l no , n £ N 



^n+n := max||am|, Re^y/ (\k + \\ - \ + n + nf + J7_ - \am\\ j, 



aIT-L := + \\ - \ + n + n f + ^+ + M- 

Proof. By Lemma 4.3, the angular operator satisfies conditions (Tl), (Bl), (B2), (Al), 
(A2), (Dl) and (D2), and the domain of the operators <Si(A) does not depend on A for 
A £ {\am\, oo). Since B* is surjective, we have 

cr ess (Si) = a ess (A) n (\am\, oo) = 
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by Corollary 2.8. Formula (44) now follows from Theorem 3.6 with C2 = \am\ and 
A e = oo. Since D is bounded and 

— \am\ \\x\\ 2 < (x, Dx) < \am\ \\x\\ 2 , 

application of Theorem 3.12 with C2 = cf = \am\ and = c\ = — \am\ yields the es- 
timates (45). By Theorem 3.6 and Proposition 3.7, no = min dimjC(_ 00) o)<S'i(A) < 

Ae(|am|,oo) 

oo. Since (\am\, X m++ i) C p(S\), the index shift no is constant in this interval. Hence 
also the assertion concerning no is proved. The explicit two-sided estimates for the 
eigenvalues X m++n are obtained if we insert the estimates (38) into (45) and observe 
that \ m+ + n > \am\ by definition. □ 

A result similar to Theorem 4.4 follows directly from standard perturbation theory 
(see, e.g., [Kat80]) applied to the angular operator with m as perturbation parameter. 
For convenience, we state this result in the next theorem. Since with the method from 
perturbation theory no index shift no occurs, a comparison of the results of the following 
theorem and of Theorem 4.4 leads to a condition for no = (see Propositions 4.8 and 
4.9). 

Theorem 4.5. Let X n be the nth eigenvalue of the angular operator A with the ordering 
described above. Then for all n £ N we have 



(\ k +\\ ~^+n) 2 + Sl_ 



/ 1 1 2 

A n < y {\k + -| - - + n) +Q + + \am\. 
The functions f2_ and 0+ are defined in (39) . 

Proof. Estimates for the eigenvalues of BB* are given in (38). Since B and B* are 
invertible, the spectrum of B = ( is given by o~ p (B) = \±^Ju^ : v n G a(BB*)}. 
Now, application of analytic perturbation theory to the operators B and A with m as 
perturbation parameter yields — \am\ < A n < y/T^+ \am\. □ 

Remark 4.6. In addition to the estimates for the eigenvalues of A presented in this 
paper, there is also another method to derive a lower bound for the modulus of the 
eigenvalues of A, see [Win06]. This method making use of sesquilinear forms yields the 
following bound: 



I An | > A Q (46) 



where 



sign(/c + \){aoo + k + |) = \auj + k + i| if au) £ [ — \k + ||, \k + \ 
2^au(k + |) if sign(|fc + \\)auj > \k + 



2i- 



In the case sign(|A; + i|)au; > \k + i| this methods yields no bound for the eigenvalues. 
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4.2 The index shift uq 



The index shift no does not depend on the choice of Ao G (|am|, A m++ i) but, of course, 
it depends on the physical parameters a, m, oo and k. 

The following lemma gives a sufficient condition for the index shift to be nontrivial. 

Lemma 4.7. If there exists an eigenvalue u of A such that 

2\am\-X m++ i < fi < A m++ i, (47) 

then we have uq > 1. If in addition A m+ +i < 3 \ am\, then there is at least one eigenvalue 
of A in [— \am\, \ am\ ]. 

Proof. Recall that A m+ +i is the first eigenvalue of A which is greater than \am\, hence 
we have u < \am\. If we also know A m+ +i < 3 \am\, then (47) shows that fi > 2\am\ — 
A m+ +i > — \am\. Hence the eigenvalue fi of A lies in [— \am\, \am\ ]. 
It remains to be shown that uq > 1. Recall that uq = min dim£/_ 00 o)Si(A) and 

\>\arn\ v ' ' 

that the right hand side is constant on the resolvent set of A and non-decreasing with 
increasing A. Hence uq = n(A) := dim£(„ oo0 )S'i(A) for all A G (\am\, X m++ i). Let V 
be an arbitrary linear manifold such that T>{S\) CDC T>(s\). In [EL04, lemma 2.5] it 
has been shown that no is equal to the dimension of every maximal subspace of 

Af{\) := {x G V : «i(A)[x] < 0} U {0}. 

Since I'(si) = D(B*) does not depend on A, we can choose T> = D(B*). Therefore 
it suffices to show that there exists x £ V(B*), x / and A £ ( | am |, A m++ i) such 
that Si(A)[x] < because then x spans a onedimensional subspace in AA(A). Since 
fj, is an eigenvalue of A, there exists an element (x, y) 1 £ V(B*) © D{B) such that 
(A-n)(y) = 0, i.e., 

{-D - n)x + By = 0, B*x + (D — fx)y = 0. 

In particular, for A G p{D) we have (D — \)~ x B*x = —(D — X)^ 1 (D — fi)y = —y + (fi — 
X)(D — X)~ 1 y and (B*x, y) = (x, By) = (x, (D + fi)x). Thus for every A > \am\ 

8i(X)[x] = (x, (-D - X)x) - (B*x, (D - X)- 1 B*x) 

= -X\\x\\ 2 2 - (x, Dx) + (B*x, y)-(fi- X)((D - X)~ 1 B*x, y) 
= (l* ~ A) {\\xg + \\y\\l) - (fi - A) 2 {(D - A)" 1 !/, y). 

Since A > \am\ = \\D\\, we haveO < - ((£>- A) l y, y) < (A-| am\) 1 ||y||2- Furthermore, 
it follows from Remark 4.1 that |x($)| = |y(vr — -d)\ for all # G (0,7r) which implies 
1 1 1 1 2 = Hz/lb- Thus we have 

Si(A)[a;] < ||x||2 (A — \am\)~ l (fi — X)(fi + X — 2\am\). 

For A := A m++ i — ^(/i + A m++ i — 2|am|) it follows from (47) that A G (|am|, A m++ i). 
Moreover, we have fi — X < and /x + A — 2\am\ = ^(fi + A m++ i) — \am\ > by (47). 
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Hence it follows that 

Si (A) [re] < ||rc||i(A- |am|) _1 (/i- A)(/x + A - 2|am|) < 0. 

□ 

Recall that v n , n G N, are the eigenvalues of BB* . 
Proposition 4.8. (i) If there exists jo > 2 siic/i i/tai 

V^o+io - V^o+io-i > 2 l am l and V^o+io+i - V^o+jo > 2|am|, (48) 
i/ien no = ?n + . 

(ii) // H-B*" 1 !! -1 > 2|am|, then the angular operator A has no eigenvalues in the 
interval [— \am\, \ am\ ] and we have no = and m + = 0. 

Proof, (i) From standard perturbation theory we know that 

sign(n)y/zv| n | — \am\ < X n < sign(n)y/i/| n | + \am\, nGZ\{0}. (49) 

Hence (48) implies that the angular operator A has exactly one eigenvalue in the interval 
[y/Vn +jo ~ \am\, ^/y no +j + \am\]. Since by (48) and (45) both \ no+ j and A m++J - lie 
in this interval, it follows that no = m + . 

(ii) Assume that ||-B* _1 || _1 > 2|am|. Then we have ^Jv\ > 2\am\ for the smallest 
eigenvalue v\ of BB*. From the estimate (49) we obtain that 

A_i < — \J~u{ + \am\ < —\am\. and Ai > ^Jv\ — \am\ > \am\ 

Hence the angular operator A has no eigenvalues in [— \am\, \am\ ] which implies m + = 0. 
For A > \am\ define the set as in the proof of Lemma 4.7. Since no is equal to the 

maximal dimension of subspaces of M{\) for A G (|am|, A m++ i), it suffices to show that 
A/"(A) = {0} for A close enough to \am\. To this end fix an arbitrary x G T^{S\) = T>(B*). 
Then it is easy to see that for all A > \am\ 

Si(A)[x] = (x, (-D - X)x) - (B*x, (D - X)~ 1 B*x) 

> (-\am\ - A) ||rc||| + (\am\ + A) -1 

> (\am\ + A)" 1 ||x||| ( H^ 1 !!" 2 - (A + \am\f). 

Since by assumption Hi?*" 1 !!" 1 > 2|am|, we have Si(A)[x] > for all x G T>(B*) \ {0} 
if A is sufficiently close to \am\. □ 

The next proposition follows immediately from Proposition 4.8. Recall that A^ +no and 
[ul 

A„_j_ no are the upper and lower bounds provided by the variational principle for the 
(m + + n)-th eigenvalue of A (see Theorem 4.4). 
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Proposition 4.9. (i) // there exists jo > 2 such that 

A !lUo - A S + io-i > ° ™* ^I + io + i - > o. (so) 

f/ien no = m+. 

(ii) 7f A^ > \am\, then the angular operator A has no eigenvalues in the interval 
[—\am\, \am\ ] and we have uq = and m + = 0. 

Note that we have 

lim \M *M 

lim A Ar+1 - A w 



= -2\am\ + lim + ±| + iV + ±) 2 + 0_ - + || + AT - ±) 2 + 0+ 

= — 2|am| + 1; 
therefore, (i) of Proposition 4.9 holds whenever \am\ < 

It can be shown that for fixed parameters a, m and u, we have lim ||-B _1 || = 0, 

|fe|— >oo 

see [Win06]. Hence, if the norm of the wave number k is large enough, then the angular 
operator has no eigenvalues in [— \am\, \am\ } and the index shift no vanishes. 



4.3 Comparison with numerical values 

Suffern, Fackerell and Cosgrove [SFC83] have obtained numerical approximations of the 
eigenvalues A of the angular operator by expanding the solution of the angular equation 
in terms of hyper geometric functions. They derived a three-term recurrence relation 
for the coefficients in the series ansatz. For the eigenvalue A they found an expansion 
with respect to a(m — uj) and a(m + to) 

A = ^C r , s a r+s {m-uj) r {m + uj) s 

r,s 

with the coefficients C T)S obtained from the recurrence relation for the eigenfunctions. 
The numerical values of [SFC83] will be denoted by A^ m ' . They differ from the eigen- 
values given in this work by a factor —1 due to the choice of the sign of the constant of 
separation in the separation ansatz for the Dirac equation. 

For fixed values of am and au>, Tables 1 and 2 contain the numerical values for the first 
positive and first negative eigenvalues A^ tabulated in [SFC83] and the analytical 

lower and upper bounds A^ and A^ from Theorem 4.4 for wave numbers k = —5, . . . , 4. 
In addition, the values for the lower bound Xq from Remark 4.6 are listed. For all 
physical parameters under consideration, apart from the case am = 0.25, auj = 0.75, 
k = — 1, we have 

WB- 1 ]]- 1 = > Re^(\k + \\ + i) 2 + fi_ ) > 2\am\ 
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where v\ is the first eigenvalue of BB* , see (38), so that we have no = and m + = 
by Proposition 4.8 (ii). Therefore, the first positive eigenvalue is indeed the analytic 
continuation of the first positive eigenvalue in the case a = 0. The case am = 0.25, 
au = 0.75, k = — 1 is discussed in the subsequent remark. 

Remark 4.10 (am = 0.25, ato = 0.75, k = -1). 

(i) In this case, the bound Xq from Remark 4.6 is not defined because of sign(fc + 
~)au) = -0.75 <-\-\k + \\. 

(ii) The inequalities (38) yield no positive upper bound for ||-B _1 || so that we cannot 
use Proposition 4.8 (ii) to conclude no = m + = 0. However, since \am\ < ^, it 
follows from Proposition 4.9 (i) that no = m + . By Theorem 4.5 we still have 

-0.25 < y/u[-\am\< X 1 <^+\am\ < 1.28078 (51) 

where Ai is the analytic continuation of the first positive eigenvalue in the case 
a = and v\ is the first eigenvalue of BB* which we have estimated according 
to (38). 

(iii) Even a positive lower bound for Ai can be obtained by means of analytic pertur- 
bation theory if a is treated as the perturbation parameter. For a = we have 
X n = sign(n)(|&; + ^| — \ + |n|) = n; hence for the given physical parameters we 
obtain 

ji-0.75 < A„ < Ji + 0.75, neZ\{0}. (52) 

In particular it follows that 0.25 < Ai. For all other values of n, however, the 
bounds An and A^ obtained from the more elaborate estimates in Theorem 4.4 
(where m plays the role of the perturbation parameter) yield tighter bounds than 
the formula above as can be seen in Table 3. 

Combining (51) and (52) we obtain 0.25 < Ai < 1.28078. 

Remark 4.11. In some cases, the bounds can be further improved. For am = 0.005 
and auj = 0.015 and k G {—5, . . . , 4} we have a(A) n [— \am\, \am\ ] = because of 
|A±i| > X^i = — \am\ > \am\. Furthermore, = -y/zTj" > \am\ so that the 

assumption of Lemma 3.38 of [WinOG] is satisfied. Hence it follows: 

(i) For k = 0, . . . , 4 we have (-(Ii? -1 !! -1 , -\am\) n a (A) = 0, hence 

A-i < H^ir 1 = — Af* 1 - \am\. 

(ii) For k = —5, . . . , —1 we have (— — \am\) n cr(A) = 0, hence 

Ai > WB- 1 ^ 1 = X [ ! ] + \am\. 
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Analogously, for am = 0.25, auo = 0.75 the upper bound for A_i can be improved if 

k = 0, . . . , 4 and the lower bound for Ai can be improved if k = —5, . . . , —2. 

Note, however, that for k = — 1 the assumptions of Lemma 3.38 in [Win06] are not 

fulfilled. 

The discussion in Remarks 4.10 and 4.11 shows that it is not easy to decide a priori 
which analytic bound gives the sharpest bound for the eigenvalues of A. It seems that 
often a combination of the various estimates yields the best result. 
It can be seen from the tables that in most cases the estimate \± yields the sharpest 
lower bound. On the other hand, Figures 2 and 3 suggest that for increasing am and 
auj the estimate Aq provides a better lower bound for the smallest positive eigenvalue 
than A^ does. 

In Figure 1 the numerical values Ag" m ' and Ag"™' together with the analytic lower 

bounds ±Aq and the analytic upper and lower bounds A^ and Aj" from Theorem 4.4 
for the lowest eigenvalues as functions of k are plotted. 



am = 0.25, auj = 0.75 





Aq 


A? 


^ [tiura] 
A S,1 


^ [ntim] 
A S,-1 


A? 


k = -5 


3.75000 


3.93330 


4.29756 


-4.34936 


4.61606 


-4 


2.75000 


2.91228 


3.30870 


-3.37371 


3.65037 


-3 


1.75000 


1.87132 


2.32657 


-2.41349 


2.71221 


-2 


0.75000 


0.75000 


1.35984 


-1.48903 


1.85078 


-1 


undefined 


(0.25000) 


0.44058 


-0.67315 


(1.28078) 





1.22474 


0.75000 


1.59764 


-1.47645 


1.85078 


1 


2.25000 


2.09521 


2.65654 


-2.57663 


2.90754 


2 


3.25000 


3.21410 


3.68229 


-3.62219 


3.93273 


3 


4.25000 


4.27769 


4.69685 


-4.64856 


4.94707 


4 


5.25000 


5.31776 


5.70622 


-5.66583 


5.95636 



Table 1: Analytic bounds and numerical approximations for the first positive and first 
negative eigenvalue of A. The estimate Aq from Remark 4.6 is a lower bound for 

|A±i|. A^ and A^' from Theorem 4.4 are upper and lower bounds for A±i. The values 
^[num] aJ^"™' for the first positive and the first negative eigenvalue of A are taken 

from [SFC83]. We have obtained the numerical values \^ um ^ and \^ m ^ D y approxi- 
mating a solution of the continued fraction equation for A. Note that for k = 0, . . . , 4 
the upper bound for A_i can be further improved, while for k = —2, . . . , —5 the lower 
bound for Ai can be improved, see Remark 4.11. For k = — 1 see the discussion in 
Remark 4.10. 
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am = 0.005, auj = 0.015 





Aq 


A? 


[niim] 
A S,1 


, [ntim] 


am 


k = -5 


4.48500 


4.97998 


4.98591 


-4.98682 


4 


99299 


-4 


3.48500 


3.97997 


3.98611 


-3.98723 


3 


99373 


-3 


2.48500 


2.97996 


2.98643 


-2.98786 


2 


99498 


-2 


1.48500 


1.97994 


1.98700 


-1.98901 


1 


99749 


-1 


0.48500 


0.97989 


0.98834 


-0.99170 


1 


00500 





0.51500 


0.99500 


1.01167 


-1.00836 


1 


01989 


1 


1.51500 


2.00249 


2.01300 


-2.01101 


2 


01994 


2 


2.51500 


3.00498 


3.01357 


-3.01215 


3 


01996 


3 


3.51500 


4.00623 


4.01389 


-4.01278 


4 


01997 


4 


4.51500 


5.00699 


5.01409 


-5.01318 


5 


01998 



Table 2: Analytic bounds and numerical approximations for the first positive and first 
negative eigenvalue of A. The estimate Xq from Remark 4.6 is a lower bound for |A±i|. 
A^ and A^ obtained in Theorem 4.4 are upper and lower bounds for A±i. The values 
^ are the first positive and the first negative eigenvalue of A calculated 
numerically by Suffern et al. [SFC83]. Note that for k = 0, . . . , 4 the upper bound for 
A_i can be further improved, while for k = —1, . . . , —5 the lower bound for Ai can be 
improved, see Remark 4.11. 
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am = 0.25 
auj = 0.75 




-5-4-3-2-10 1 2 3 4 

wave number k 



am = 0.005 
auj = 0.015 




-5-4-3-2-10 1 2 3 4 

wave number k 



A s,i ' A sw ±A i ' ±A i ±A Q 

Figure 1: The plots show the lower bound Xq for the absolute value of the eigenvalues 

of „4 and the analytic upper and lower bounds and from Theorem 4.4 for two 
different values of (am, auS). In addition, the numerical values for the first positive 
and first negative eigenvalue of A from [SFC83] are plotted. The analytic bounds have 
not been plotted in the interval (—1,0) because for wave numbers k in that interval 
the angular operator is not uniquely defined as a selfadjoint operator. Note that for 
(am,auj) = (0.25, 0.75) the bound Aq is not defined for k £ (—1.25, 0). In the case 
(am,auj) = (0.005,0.015) the analytic lower and upper bounds A^ and Aj"' are so close 
to each other that they seem to coincide in this resolution. 
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am = 0.25, auj = 0.75 



7, 

A 


n 

— U 


An 


, [num] 
A S,n 


x [num] 
A S-n 


An 


n 


= 1 


0.75000 


-1 f f~\ r—rr-' A 

1.59764 


-1.47645 


1.85078 




2 


1 75000 


2 22587 


-2 93549 


2 60850 




3 


2 75000 


3.17408 


-3 16265 


3 50000 




A 

-k 


O. f ouuu 


A 1 ^1 97 


A 1 9AA(\ 


A AA(\7(\ 




K 
■J 


A 7^000 


o. LyJoOO 


0. 1UUOO 








x[J] 
An 








h — 


1 

— ± 


, [num] 
A S,n 


x [num] 
A S-n 


An 


n 


= 1 


(0.25000) 


0.44058 


-0.67315 


(1.28078) 




2 


1.33114 


1.84225 


-1.87948 


2.26556 




3 


2.48861 


2.90717 


-2.92301 


3.26040 




4 


3.55789 


3.93475 


-3.94336 


4.25780 




5 


4.59768 


4.94973 


-4.95513 


5.25625 



Table 3: For am = 0.25, aw = 0.75 and k = 0, —1 the numerical values Ag_ n and 

Aj™' and the lower and upper bounds An and A^ from Theorem 4.4 are shown. For 
k = —1, n = 1 we refer to Remark 4.10. 



am = 0.005, au = 0.015 



k 


= 


\M 
An 


x [num] 
A S,n 


x [num] 
A S-n 


X M 

An 


n 


= 1 


0.99500 


1.01167 


-1.00836 


1.01989 




2 


1.99500 


2.00435 


-2.00369 


2.01249 




3 


2.99500 


3.00273 


-3.00245 


3.01000 




4 


3.99500 


4.00180 


-4.00184 


4.00875 




5 


4.99500 


5.00158 


-5.00148 


5.00800 






An 








k = 


-1 


< [num] 
A S,n 


x [num] 
A S-n 


X M 

An 


n 


= 1 


0.97989 


0.98834 


-0.99170 


1.00500 




2 


1.98749 


1.99567 


-1.99636 


2.00500 




3 


2.99000 


2.99730 


-2.99759 


3.00500 




4 


3.99125 


3.99803 


-3.99819 


4.00500 




5 


4.99200 


4.99845 


-4.99855 


5.00500 



Table 4: For am = 0.015, aw = 0.025 and k = 0, —1 the numerical values A^_ n and 
An™' and the lower and upper bounds An and A^' from Theorem 4.4 are shown. 
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k = 
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m = 


0.025 
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0.75 
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-2.5 



■1.5 -1 -0.5 0.5 1 



1.5 



k 
m 
u 



4 

0.025 
0.75 



■1 



Affl 



Figure 2: Lower bounds for the eigenvalue of A with smallest modulus for k, m and 
uj fixed. Note that for both cases of k, for each of the plotted estimates there is an 
interval for a where it provides a larger lower bound for the modulus of the eigenvalues 
of A than the other bound. 
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Figure 3: Lower bounds for the modulus of the eigenvalues A of A for m = 0.25 and 
u> = 0.75 fixed. In the first graph, the bounds are plotted as functions of a for k = 
fixed. The second graph shows the bounds as functions of k with a = 1 fixed. Note 
that physically only the values for integral values of k make sense. 
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